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Abstract 

Vacuum polarization in QED in a background gravitational field induces interactions which 
effectively modify the classical picture of light rays, as the null geodesies of spacetime. These 
interactions violate the strong equivalence principle and affect the propagation of light leading to 
superluminal photon velocities. Taking into account the QED vacuum polarization, we study the 
propagation of a bundle of rays in a background gravitational field. To do so we consider the 
perturbative deformation of Raychaudhuri equation through the influence of vacuum polarization 
on photon propagation. We analyze the contribution of the above interactions to the optical scalars 
namely, shear, vorticity and expansion using the Newman-Penrose formalism. 
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I. INTRODUCTION 



Vacuum polarization is an essential ingredient of QED whose contribution leads to aston- 
ishingly precise agreement between predicted and observed values of the electron magnetic 
moment and Lamb shift. On the other hand being a quantum field theortetic effect it would 
be interesting to look for its implication in semi-classical gravity as a quantum characteristic 
of the electromagnetic field coupled to the gravitational field [l| ■ Effects of QED interactions 
such as vacuum polarization on the photon propagation in a classical background gravita- 
tional field gives rise to a wide range of new phenomena. Simple analysis of null rays in 
classical general relativity implies that a curved spacetime, compared to the flat case, could 
be treated as an optical medium with a refractive index [2,]. Now adding to that QED effects, 
in the context of semi-classical gravity, leads to interesting phenomena such as dispersive 
effects and polarization dependent propagation of photons. One of the main consequences 
is found to be the light cone modification in such a way that the QED photons do not prop- 
agate along the null geodesies of the background geometry. They propagate instead along 
null geodesies of an effective geometry. In many cases depending on the direction and po- 
larization of photons, superluminal propagation becomes possible. Quantum characteristics 
of photon propagation in a curved background opens up a whole new field coined quan- 
tum gravitational optics (QGO)[3], Started with the pioneering paper by Drummond and 
Hathrell who considered the effect of vacuum polarization on photon propagation in a curved 
background and continued with detailed study of the same effect for specific spacetimes {4 1- 
s|. Coupling of the electromagnetic field with curvature in the QED action introduced in 1] 
violates the strong equivalence principle (SEP) in a mass scale comparable to the electron 
mass m. The modification of the underlying spacetime geometry felt by photons due to QED 
interactions has some unexpected results. The Newman-Penrose (NP) formalism has been 
found to be an elegant way of analyzing the results in the SEP violating cases. Quantum 
gravitational optics attributes velocity shift and null cone modification at each point to only 
a single NP scalar for each of the Ricci and Weyl tensors, namely $00 and \l/o respectively 
0]. The question is what would be the role of other NP scalars in determining the forms of 
Ricci and Weyl tensors? or what, if any, would be the effect of the gravitational field on the 
other photons traveling with the normal speed of light Here we shall study the geometric 
properties of physical null congruences which are families of physical null rays determined 
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by the effective metric G^ v . These are then distinct from the propertries of the geometric 
null congruences characterized by the spacetime metric g^ v . These congruences are special 
in the sense that the modification of the light cone in a local inertial frame (LIF) can be 
described as (?7( a )(6) + « 2cr (a)(6) (R)) k^k^ = 0, where i]{a)(b) is the Minkowski metric, a is 
the fine structure constant and C(o)(6) (R) depends on the Riemann curvature at the origin 
of the LIF. In one-loop approximation, these photons travel with unit velocity. For on-shell 
photons, some of the NP scalars play the role of optical scalars and their propagation, along 
the rays, is governed by the known Raychaudhuri equation. We show that these scalars must 
be effectively modified to describe the physical null congruences and find that a new set of 
NP scalars contribute in the definition of these effective optical scalars. The modified Ray- 
chaudhuri equation governing their propagation will be discussed. The effective geometry 
may have features that are seen only by photons when propagating in the geometry. There 
may be hidden singularities in the form of geodesic incompleteness which can be studied 
through the implementation of the effective Raychaudhuri equation. 



II. QUANTUM GRAVITATIONAL OPTICS 



The effect of one-loop vacuum polarization on the photon propagation in a fixed curved 
background spacetime is represented by the following effective action derived by Drummond 



and Hathrell 
r = / olx 



1] 



-g 



-F^F^ + [aRF tiV F" v + bR, u F» x F\ + cR^F** F") 



Here, a = — ttt-, b = M^- and c = — where a is the fine structure constant and 
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m is the electron mass. The notable feature in the above action is the direct coupling of 
the electromagnetic field to the curvature tensor which in effect violates the strong equiv- 
alence principle. When we consider the equations of motion derived from this action the 
Bianchi identity does not change but it gives rise to curvature dependent modifications to 
the equation of motion in the form 
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D,F^ + — (2bR» X D^F X » + AcR^D^o) = 0. 



(2) 



There are some approximations under which this equation of motion was obtained. The first 
one is the low frequency approximation in the sense that the derivation is only applicable to 
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wavelengths A > A c . By this approximation we ignore terms in the effective action involving 
higher order field derivatives. The second is a weak field approximation for gravity. This 
implicitly means that the wavelengths X c « L are considered, where L is a typical curvature 
scale. 

The characteristics of the light propagation can be studied by applying the geometric optics 
approximation. The electromagnetic field is written in the form = Aa^e 10 , where A 
is the amplitude and a M is the polarization vector. The amplitude is taken to be slowly- 
varying in comparision with 9. In each small region of space, we can speak of a direction of 
propagation, normal to a surface at all of whose points the phase of the wave is constant. 
If we now identify the wave vector as k^ = 8^9, the equation governing the components of 
this vector field for a photon with spacelike, normalized polarization vector i.e, a^a^ = — 1, 
can be read from (j2J) as 

k 2 - ^-R^k x + ^R^ Xp Wk x a v a p = 0. (3) 

Eq. ([3]) is an effective light cone equation, re-expressed in terms of the Weyl tensor is given 
by 

k 2 - {2b \ 4c) R, x k»k x + ^C^#»fcVa' = 0. (4) 

The corresponding momentum of photon, p**, is the tangent vector to the light ray x M i.e, 
p^ = (s). Equation (J3J being quadratic and homogenous could be written in the 
following form 

G^k,k u = (5) 

This equation implies a nontrivial relation between k^ and through Q by the following 
argument. Refer ing to the defitions of k^ and p^, the former defined as the gradient of a 
phase, constitutes the component of a one-form and hence belongs to a cotangent space 
whereas the latter being a tangent vector belongs to the tangent space. In the usual free 
electromagnetic theory they are related through the metric tensor i.e k^ = g iiv p v . In the 
modified theory however, using equations (j3J) and (jSJ), we arrive at the following non-trivial 
relation, 

K = G^tf = ?v + ^ [(6 + 2c) R^p x - (4c) C^/aV] (6) 

where G = Q^ 1 . This shows that at this level of approximation, QGO can be characterized 
as a bimetric theory. The physical light cones are determined by the effective metric G^, 
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and are distinct from the geomerical null cones which are fixed by the spacetime metric g^ u . 
(Indices are always raised and lowered using the spacetime metric g^)- 
It should be noted that all these equations are manifestly local Lorentz invariant though 
the presence of the explicit curvature coupling in the effective action means different dy- 
namics in the LIFs at different points in spacetime. In this sense these equations violate 
the strong principle of equivalence (SEP). Some of their implications including gravitational 
birefringence and superluminal speed of light have been discussed in [3] . Equation (jSJ) can 
be employed to examine the SEP violating effects which arise in the physical photon con- 
gruence. In the present study we will focus on the modifications produced by the vacuum 
polarization effects on the propagation of light rays in a curved background. 



III. EFFECTIVE OPTICAL SCALARS 



The physical congruence is specified by the integral curves 7 of the vector field k^ 
parametrized by the parameter u and scaled such that V^w = 1. In order to examine 
the relation between the neighboring curves of the congruence, the connecting vector field 
g M is introduced 8j. Defined along a particular curve 7, it characterizes the displacement 
from a point P E 7 to a point P on a neighboring curve, where P and P have the same 
parameter value u. Mathematically, this means that q^ 1 is being Lie transported along the 
curve by the vector field k^ . 

C k q^ = that is V k q M = V q k^. (7) 

Here the notation Vx = X a V a for directional covariant derivative is used. In a SEP 
violating theory like QGO, it is illuminating to work in the NP tetrad formalism. Following 
closely the notations used in [7|] and [8|, the first step is to choose a null tetrad. We choose P 
be a null vector along the photon momentum. Let a M and be spacelike transverse vectors 
(in the present case they will be identified with the polarization vectors later) and define the 
complex null vectors m M and by = {a? + ib^ 1 ) and mf 1 = {a? — iU 1 ). Adding to 
this set another real null vector n M , they satisfy the usual NP orthogonality conditions 

l.m = l.m = n.m = n.m = 0, (8) 

and 

1. 1 = n.n = m.m = m.m = 0, (9) 



m 



e (4 ) 



The normalization conditions 

l.n = —m.m = 1 

are also imposed. We assign a vierbein e M , * to this tetrad as follows [1 

e(i) = e (2) = I e(2) = e (1) = n e( 3 ) = 
with the frame metric of the form 

/ 1 \ 
10 
-1 
\0 0-l J 



-e (3) = m, 



(10) 



:n) 



(12) 



Weyl tensor being trace free is completely specified in this basis, by the following five complex 
scalars Q], 



*0 = 


-C abcd l a m b l c m d 


= — Cl313 


*1 = 


-C abcd l a n b l c m d - 


= — Cl213 


*2 = 


-C abcd l a m b nfn d 


— — Cl342 


*3 = 


-C abcd l a n b m c n d 


= — Cl242 


^4 = 


-C abcd n a m b n Ci m 


i r~i 
— —^2424 



(13) 

The ten components of the Ricci tensor are also defined in terms of the following four real 
and three complex scalars. 



$00 = 


-\Ru\ 


$22 = 


-^22; 


$02 




— R11 
2 


$20 = 


-\« 


$01 = 




$10 = 


~Ru; 


$12 




~2 R23 


$21 = 




$11 = 


-\ ( R » 


+ R34) ; 




R = 


1 

12 


(R12 - 


-R34) • 





44) 



21: 



(14) 

In tetrad formalism, various quantities called "Ricci rotation coefficients", 7( c )( a )(&), are em- 
ployed to account for the tetrad covariant differentiation. These quantities can be defined 

as 

e(a)w = e { %~f {c)ia )( b )e {b l. (15) 
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In NP formalism, they are called "spin coefficients" and are designated by the following 
special symbols [9|: 

2 (T(2)(i)(i) +7(3)(4)(i)) ; 

(7(2)(1)(2) + 7(3)(4)(2)) J 



K 


= 7(3)(i)(i); 


p 


= 7(3)(i)(4); 


e 


a 


= 7(3)(i)(3); 


(I 


= 7(2) (4) (3); 


7 


A 


= 7(2) (4) (4); 


T 


= 7(3)(i)(2j; 


a 


V 


= 7(2) (4) (2); 


7T 


= 7(2)(4)(i); 


(3 



2 
1 

2 
I 

2 



(7(2) (1) (4) +7(3) (4) (4)) ; 



(16) 



A general rule to obtain the complex conjugate of any quantity is to replace the index (3), 
wherever it occurs with the index (4) and vice versa. 

Using the NP formalism has some advantages, including the fact that entirely scalar quan- 
tities are used. Also a considerable economy of notation is achieved when complex spin 
coefficients are used instead of the Christoffel symbols employed in the conventional coor- 
dinate approach. This formalism is most advantageous if the null tetrads introduced could 
be completely tied to the geometry of the problem. Otherwise some freedom remains in 
the choice of the basis as a result of subjecting it to the Lorentz transformation. This has 
the effect that many of the quantities involved in the calculation are of the same nature 
as gauge quantities whose values transform in certain ways as the basis frame is varied in 
accordance with the remaining freedom. In QGO, one null direction is singled out and that 
is the unperturbed photon momentum which is fixed in the direction of the null vector P. 
Therefore the allowed transformation would be those which leave the P direction unchanged 
and preserve the underlying orthogonality and normalization conditions. These are classified 
as follows: 



class / : / i — > I; m i — > m + al\ m i — ► m + a*l; n i — > n + a*m + am + aa*l; 
class II : I i — > A~ l l; n i — > An; m i — > e l9 m; m i — > e~ l9 rn, (17) 



where a is a complex function and A and 9 are two real functions on the manifold. We shall 
now proceed to write down equation (j7j) in the above tetrad basis. Before doing so we employ 
all the gauge freedom we have to fix the tetrad basis. As the tetrad frame is transformed, 
different tetrad components are subject to changes. With / vectors as the velocity vector 
along the geometric congruence, we have 



k = 0, 



(18) 



and if they are afnnely parametrized 

V,P = k = e = 0. (19) 

If k = 0, the latter requirement can be met by a class II rotation which will not affect the 
direction of I nor of initially vanishing of k. By a suitable rotation of class I, it could be 
arranged so that tt = 0. So the gauge we are working in is the one in which k = e = n = 0. 
After such a rotation, the newly oriented vectors n, m and m will remain unchanged as they 
are parallely propagated along I. This could be invoked through the following relations 

S/ im u = ttT - kti v + (e - e*) m u 

V,n" = ir*m u + 7im u - (e + e*) n u (20) 

which vanish in the above gauge. Now we have exploited all the gauge freedom we had and 
ready to write equation ([7]) in the above gauge. We only consider the states which propagate 
with a well defined polarization 181 ] . For two transverse polarizations expressed in terms of 
the vectors and m^, representing the left and the right handed circular polarizations, we 
have 



After multiplying by e£ , we get 



l v + \ (6 + 2c) R Xu l x =F ^C\ CR l u (m x ± m A ) {m K ± m K ) 



(21) 



Vkq (c) = j v (c)(a) + -La^M ) q {b h(a K l Kb ) + g ( a)7 (a)(C)( %), (22) 



where 

A (c)(a) ;= + 2c) R {c)(a) =f (2c)M (c)(a) , 

M (c)(a) ._ C f(c)(3)(a)(3) + (c) (4) (a) (4) ± ^ (c)(3) (a) (4) + ^nr (c)(4) (a)(3) j _ (23) 

The specifically gravitational birefringence shows up here in M^ a ^ c \ (we note that = 
A(°)( c )). Equation ([22]) is a system of coupled differential equations describing the propaga- 
tion of the tetrad components of the connecting vector along the velocity vector k^ 1 . The 
second term in f[2"2"]) appears as a consequence of the null tetrad propagation along the phys- 
ical ray (with the tangent vector A;^) i.e q u V k^f 1 ■ We introduce different tetrad components 
of q^" in the following way 

(f = g F + + lm» + hn". (24) 



so that 

q (1) = q(2)=g; q {2) = q W = h; = -q {A) = £; = -q (3) = £. (25) 
Explicitly for c = 4, ( |22i) gives 

= + i,^) [g^ 7{ 3 )(1){2) + q (3 hmm + q ( %mm] 

+ i A(4K4) b (2) 7(4)(l)(2) + g (3) 7(4)(l)(3) + 9 {4) 7(4)(1)(4)] 

+ i^VWw - ^ 2){ % {a h^m (26) 

/ / i lib 

in which the last term is again the contribution due to the null tetrad propagation along 
the physical ray discussed below equation (1231) . To identify the geometrical and physical 
content of this equation, a "covariant approach" in which any complicated gauge behavior is 
completely avoided must be chosen. For example, in (I26I) . the term A^^q^j^^^) having 
such a complicated gauge behavior, transforms into 

[ A WW - 2aA^ + a 2 A^} q^ [ 7(4)(1)(2) + a* 7 (4)(i)(3) + a 7 (4)(i)(4)] (27) 

under a class / transformation. To keep the formalism covariant we remove the gauge- 
dependent terms by applying the extra condition that the neighboring pair of rays to satisfy 
the following condition, 

q.k = 0, (28) 
and call them abreast . This is a primary constraint and it means 

O (a ) : q.l = =>- h = 0; 

O (a 1 ) : A ( i )(o/2) = $oo = $oi = $io = = #i = 0. (29) 

Eq. (I2"9"j) expresses the abreastness constraint which is independent of the parametrization 
of k^ as could be easily checked by making the substitution g M i— ► q^ + Ak^. It is easy 
to see that under the new parametrization, q.k = 0(a 2 ), and therefore could be ignored. 
Furthermore, like any other constrained system, we need to make sure that the eq. (|28|) is 
conserved under the propagation along k. In other words, (q-k) = does not lead to a 
secondary constraint. Up to the first order in a, we have 

V k = =>- (V 9 P) k„ + <f (V fc Ag = 

O (a ) : (V g P) Z„ + g" (W M ) = 

O (a 1 ) : [A«< a ty 6) + A«<V a) ] 7(o)(1){6) = 0, (30) 
9 



which leads to the same constraint mentioned above. We emphasize that provided the 
abreastness constraint fl29l) is satisfied, the tangent vector k^ can be described as 



k(2) = 1 J7v4 (1 )( 2 ), &(a^2) = 0, 



rn 



k 2 = (a 2 ) , V fc fc" = O (a 2 ) 



(31) 



The abreast rays are special physical rays whose effective light cone condition, in a local 
inertial frame, is given by 



(V{a)(b)+a 2 a ia}ib) (i2))*W 6) = 0. 



(32) 



We recall from ([3]) that the general light cone modification induced by the effective action 
([1]) is of order a. 

For abreast rays, the covariant deformation of £ implies 



(33) 



?7r rrr 

Due to their gauge-dependence and also for our purpose, propagation of the other compo- 
nents of the connecting vector q are of less geometric and physical importance and we just 
state the results below for the sake of the completeness of the discussion, 



V k h = 



VkQ = £ (P + a*) + £ (a + (3*) + -L [A (1)(3) (£p + £<r*) + c.c] 



(34) 



Now having the equation (1331 we are ready to study the evolution of the cross sections of 
(physical) bundles of rays, formed by their intersection with the spatial plane spanned by 
m M and m M . To do so we employ the so called optical scalars associated with a bundle of 
rays, namely the expansion, the shear and the twist 10j. Since the transverse distances 
from a specified ray 7 are (frame)observer-independent, for abreast rays these scalars bear 
explicit physical interpretations. To simplify the case we consider the propagation along k 
of the area of a small triangle formed by the points o, £1 and £2 (see figure 1), we have 



V k SA = V fc 
The scalar quantity 



(£i£ 2 - 6£i) 



p 



-L(A^p + A^a)+c.c 
m 2 



5A 



-Re 



p --L(A^ P + A^a) 



(35) 



(36) 
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mm-plane 



FIG. 1: Three abreast rays of a congruence and the small triangle formed by two connecting vectors 
£1 and £2; connecting the ray 7 to the rays 71 and 72 in the mm-plane. 



is called the expansion parameter, 6 e g, whose role as a measure of the pattern convergence 
(or divergence) is clear from (1351) . This equation shows that in an effective theory, unlike in 
the classical case, the local effect of a can also change the area. Another useful parameter, 
the "luminosity parameter" L, is defined along a bundle of rays as L 2 oc 5 A and is related 
to 9 e s as follows (using (1331)). 

V k L = 6 eS L. (37) 

In other words, the expansion parameter is the logarithmic derivative of the luminosity 
parameter. An alternative formula for (9 e ff is given by ifc** 
Setting 9 ef[ and the coefficient of £ equal to zero in fl33l . we obtain 



V fc £ 



-ilm 



p - -L (^(3)(4) +A (3)(3) a \ 

m 2 



(38) 
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So it could also be claimed that Im [p — ^ (A^^p + A^ 3 >^ 3 'aj\ measures the twist in the 
bundle's cross section and could be called the effective twist, u e s . The fact that this combi- 
nation of scalars is a measure of twist is consistent with the requirement that it measures the 
failure of fc M to be hypersurface orthogonal. In other words it should vanish for a congruence 
corresponding to the gradient of a scalar field i.e, 

K = V M / =4> u cS = (39) 

To show this, we note that the left hand side of (|39|) is equivalent to Vdj = 0. At the 
leading order, 0(a°), this gives 7( a )(i)(&) = 7(&)(i)(a), which in turn, using ( fl6l) . leads to 
Imp = for a = 3 and 6 = 4. The subleading term, O (a 1 ), is given by 

[(& + 2c) R Xu - (4c) C vpXa a p a°\ e {a) \ a){l){b) ef = [(& + 2c) R XfX - (4c) C wA(T </a CT ] e (e)A 7(e ) {1)(/) e 

(40) 

Multiplying both sides of (00]) by e ^ d) e y ^ c) e p{g) e a{h \ we end up with 

A (a)(c) 7(a)(i)(^ (6)(d) = ^ (e)(d) 7(e)(i)(/)^ (/)(c) (41) 

Using the abreastness conditions (|29|) . for c = 4 and d = 3, we obtain 

, proving our assretion that u c g = to the right order. 
If # e ff = uj ef[ = 0, the remaining part of eq. (|33|) . i.e 

1 



mr 



(AWa + AWp*), (42) 



can be interpreted as the effective shear, cr e fj, which is a measure of the distortion in the 
shape of the bundle's cross section so that a circular cross section transforms into an elliptic 
one. Setting a c s = se 2W , ( 1331 reads Vfc£ = — se 2i9 ^. This shows that, Vfc£ is a real multiple 
of £ when arg£ = 9, 9 + % or 9 ± |7r, where for s > 0, in the first two cases we get contraction 
towards the origin while in the other two cases it experiences dilation. Thus s is a measure 
of the pattern's shear while 9 defines the angle that its minor axis makes with the £-plane. 
The quantities u e s and a e s can also be defined as 

kc ff | 2 + fl e 2 ff = \k^. (43) 
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IV. PHYSICAL NULL CONGRUENCE AND SPACE TIME CURVATURE: EF- 
FECTIVE RAYCHAUDHURI EQUATION 



Since the evolution of the cross section of a bundle of rays is characterized by the quantities 
p and a in the one hand and on the other hand the spacetime curvature affects the geometry 
of null congruences, we need to study the variation of optical scalars along the physical ray 
. This can be achieved by examining the second derivative of the connecting vector, q^, 
propagating along the wave vector through the operation of on eq. ([T]) , 

V fc V fc ^ = V fc (V,F) = V,V*fc" + R Xu fk x q»k° = R Xu /k x q»k° + O (a 2 ) . (44) 

After multiplying by and taking the components corresponding to the parallely propa- 
gated tetrads we get 

V k V k qU = i? (1)(3)(1) (c) e + ^ (1 )(4)(i) (C) e + O (a 2 ) . (45) 

Here, we have neglected the terms including the product of curvature components which 
would be suppressed by O(^). From ( 1331) . we obtain the useful (2 x 2) matrix form as 

V fc Z = PZ, (46) 

where 



For abreast rays (130]) . the equation (T44]) gives 



(47) 



V fe V fc Z = 0. (48) 
After differentiating fj46|) once more and using fT48l . we get 

V fc V fc Z = (VjfeP) Z + P (VjfeZ) = (V fe P + P 2 ) Z = 0. (49) 
Since this holds for arbitrary £, we must have 

V fc P = -P 2 . (50) 

Written out in full, eq. (1501) gives the directional propagation of p and a along the physical 
ray as follows, 
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Vfc p = p 2 + \a\ 2 + JL (p 2 + |a| 2 ) + A^pa + A^ pa*] 

VkG = (p + p*) a + -L [A^W (p + p*) a + A<™ V + A (4)(4) |p| 2 ] (51) 

Compared with the results in the classical case, the above equations could be called effective 
Sachs equations. With these equations in hand, we can derive the equations governing the 
directional propagation of the effective optical scalars as 

VfcCT cfr = — 20 e ffCTeff, (52) 

Vfc^eff = — 29 cS uj e ff, (53) 
V k 6 cS = u 2 eS - e 2 ^ - kff| 2 • (54) 

These equations are the main results of the paper and their physical importance (conse- 
quences) is discussed below. 

V. DISCUSSION 

Here we employed the Newmann-Penrose tetrad formalism to derive the effective Ray- 
chaudhuri equation, i.e the Raychaudhuri equation modified by the vacuum polarization 
effects on the propagation of a bundle of rays. In the derivation of the classical version of 
the Raychadhuri equation [ll|-[l3] the variations of the tangent vector k^ v is projected [3] 
on the plane spanned by m fl and m M and then its trace , antisymmetric and symmetric parts 
are designated as expansion, vorticity and shear respectively. Here in the modified version 
we chose, through the equation (J7]), an equivalent method in which the variations of the 
connecting vector q 11 instead of the tangent vector have been considered and then defined 
the corresponding effective optical scalars . The invariants of the theory are easily seen in 
this way. For a luminosity parameter, L, we see from (1531) that (L 2 u e ff) = 0, in other 
words here q.k and L 2 uj e ff are the constant geometric quantities along the congruence. For 
two neighbouring rays of 7, whose connecting vectors q^ and q^ independently satisfy (J44j), 
a symplectic invariant is attributed as follows 

g M V fc ^ - q»V k q„. (55) 

These invariants reduce to the classical ones in the limit of zero perturbation and therefore 
there remains no anomaly in QGO. This is the case since we have applied the symmetries 
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present at O (a ) as a set of constraints and exploited them to find the covariant quantum 
corrections. Since the effective action we started with has been derived in a gauge invariant 
manner the results obtained are also gauge invaraiant. 

The classical limit of the eq (j53|) can be compared with that corresponding to the standard 
(pure general relativistic) Raychaudhuri equation, namely 

V fe # = u 2 - 9 2 - H 2 - $oo- (56) 

The last term, i.e. $00 — f^WC 1 )' * s ne g a tive for all known forms of ponderable matter. 
Then for an initially contracting congruence, in the absence of rotation and shear (p real 
and a = 0), the eq. fl56l) results in the divergence of the expansion parameter (9 — > —00). 



This is a necessary though not a sufficient condition in singularity theorems 13]. In our case 
a real p and vanishing a is equivalent to u e s = a c s = 0. Using (ED), we see that during the 
propagation along the ray, they differ from zero only by O (a 2 ), 

V k d = Vjfcfleff - 0(3)(4)P) 

= + ^A3)(4)# 2 + O (a 2 ) (57) 

On the other hand, the last term in ( 156]) does not appear in our calculations due to the 
abreastness conditions. However, even if Rmii) = 0, the correction terms calculated in this 
paper will contribute and their sign must be taken into account. In QGO, the light cone 
condition as well as the velocity shift depend on $00 and \I/o ( or -4(i)(i) in our notation), 
while in the propagation of the expansion parameter these are $n, A and ^2 (or ^4(3)(4) in 
our notation) that play role purturbatively. These may have consequences including exis- 
tence of singularities in the effective geometry that are seen only by QGO photons. The 
existence of singularities seen by nonlinear photons has been studied in [141 ] . The gener- 
alization of Drummond-Hathrell results to massless neutrino propagation in a background 
spacetime 15] can be employed to find the corresponding effective Raychaudhuri equation. 
On the other hand, the effects originated from QED interactions discussed in this paper, 
in which local invariance and weak equivalence principle are preserved, are somewhat simi- 
lar to the phenomenological theories incorporating Lorentz and CPT violating interactions. 



For example, in a theory describing the propagation of light discussed in 16|, the Lorentz 
violating coupling constant K^ v \ p plays the role of the Riemann tensor R^ v \ p in eq. ([T]). 
It is therefore easy to translate the phenomenological aspects of the effective Raychaudhuri 

15 



equation discussed here to this class of theories. We recall that the relative shift in velocity 
or optical scalars would be of O for a typical curvature scale L. 
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